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ABSTRACT 

Three observations of the 5.54 s Transient Anomalous X-ray Pulsar XTE J1810— 197 obtained over 6 
months with the Newton X-Ray Multi-Mirror Mission (XMM-Newton) are used to study its spectrum 
and pulsed light curve as the source fades from outburst. The decay is consistent with an exponential 
of time constant w 300 days, but not a power law as predicted in some models of sudden deep crustal 
heating events. All spectra are well fitted by a blackbody plus a steep power law, a problematic model 
that is commonly fitted to anomalous X-ray pulsars (AXPs). A two-temperature blackbody fit is also 
acceptable, and better motivated physically in view of the faint optical/IR fluxes, the X-ray pulse 
shapes that weakly depend on energy in XTE J1810— 197, and the inferred emitting areas that are 
less than or equal to the surface area of a neutron star. The fitted temperatures remained the same 
while the flux declined by 46%, which can be interpreted as a decrease in area of the emitting regions. 
The pulsar continues to spin down, albeit at a reduced rate of (5.1 ± 1.6) x 10~ 12 s s _1 . The inferred 
characteristic age t c = P/2P w 17,000 yr, magnetic field strength B s » 1.7 x 10 14 G, and outburst 
properties are consistent with both the outburst and quiescent X-ray luminosities being powered by 
magnetic field decay, i.e., XTE J1810— 197 is a magnetar. 

Subject headings: pulsars: general — stars: individual (XTE J1810-197) — stars: neutron — X-rays: 
stars 



1. INTRODUCTION 

The bright 5.54 s X-ray pu lsar XTE J181 -197 was 
discovered serendipitously by llbrahim et alJ <|2004t) us- 
ing the Rossi X-ray Timing Explorer (RXTE), and was 
localized in two target-of-opportunity (TOO) observa- 
tions with the Chand ra X-ray Observatory, reported by 
IGotthelf et alJ 1)2004 hereafter Paper I) and llsrael et all 
( 2004) . The source had became active sometime between 
2002 November 17 and 2003 January 23. The maximum 
(2-10 keV) flux observed by RXTE was w 6 x 10" 11 
ergs cm -2 s _1 , when it was already declini ng with an ex- 
pone ntial time constant of 269 ± 25 days ijlbrahim et al.l 
2004). Of unique value are all of the archival Ein- 
stein, ROSAT, and ASCA detections at the location 
of XTE J1810— 197, which indicate a long-lived quies 
cent baseline flux of F x (0.5 - lOkeV) » 7 x 10 
ergs cm~ 2 s _1 for at least 13 years and possibly for 
23 years prior to the onset of the active state (Paper I 
and Figure ^) . Timing of the highly modulated signal 
over the first 9 months of the outburst indicated rapid 
spin-down with characteristic age r c = P/2P < 7600 yr, 
surface magnetic field B s w 2.6 x 10 14 G, and spin-down 
power £k 4 x 1 33 erg s" 1 l lbrahim et al.ll2004[) . Deep 



13 



2004) detected a faint can- 



IR observations ijlsrael et al.l 
didate of K s = 20.8, similar to ones associated with other 
AXPs, within the final 0."6 radius Chandra error circle 
of XTE J1810-197. Variability of the IR s ource con- 
firmed it s identification with XTE J1810-197 l|Rea et alJ 

2004©. 

The magnetar model ; Du ncan fc Thompsonl Il992j) 
seems poised to address observational properties of sev- 
eral categories of single, young neutron stars (NSs) that 
are not powered by rotation. These include the AXPs, 
Soft Gamma-ray Repeaters (SGRs), and Central Com- 



pact Objects (CCOs) within supernova remnants that 
lack detected periods but whose spectra resemble those 
of AXPs and SGRs. Such unification is supported 
by the discovery of transient behavior in AXP-like ob- 
jects. The first suggested TAXP was AX J1844.8-0256 
ijGotthelf fc Vasishij Il998|l . The timing and spectral 
properties of that 7-s X-ray pulsar are consistent with 
an AXP, while its flux was found to change by orders- 
of-magnitude from an "active" pulsar state to a faint 
"quiescent" state since its discovery. However, the lack 
of a P measurement for AX J1844.8— 0256, which was 
observed serendipitously only once, has left its status 
as an AXP unconfirmed. With the discovery of the 
bona fide TAXP XTE J1810-197, consideration of pos- 
sible related observations such as th e cooling trend in 
the p ost-burst fluxes of SGR 1627-41 llK o uyeliotou et alJ 
Hill, SGR 1900+14 llLenters et all 12003ft . and AXP 
IE 2259+586 <|Woods et al.l2004ft. and the variable CCO 
in th e young SN R RCW 10 3 iGotthelf .Petre. fc Hwand 
119971 IGotthelf. Petre. fc Vasishll Il999aft . motivates de- 
velopment of a unified magnetar model for the entire 
AXP/SGR/CCO family. 

Although TAXPs are relatively rare, their short ac- 
tive duty cycle suggests the existence of a significant 
number of as-yet unrecognized, although not necessar- 
ily undetected, young NSs. XTE J1810-197 provides 
a new window into this population, with a well docu- 
mented quiescent history and detailed observations dur- 
ing its active, pulsed phase. Monitoring of the fading 
of XTE J1810— 197 takes advantage of a possibly unique 
opportunity to study evolving spectral and timing behav- 
ior of a TAXP, which reveals the time-dependence of the 
thermal and non-thermal emission mechanisms to which 
magnetars convert their energy. In particular, the physi- 
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Fig. 1.— The 0.5 - 10 keV flux history of XTE J1810-197 
spanning 24 years. The historical data points are from Paper I. 
The solid line is an exponent ial fit with time con stant 269 days 
to the RXTE monitoring from Ibrahim ct al. (20041). Upper limits 
indicated in 2002 are also from several months of monitoring by 
RXTE, which constrains the onset of the outburst to < 2 months 
before its initial detection in 2003 January. 



cal interpretation of the two-component X-ray spectrum 
that seems characteristic of most AXPs can be investi- 
gated in more ways than is possible for a quiescent source. 
We are beginning such a study with the XMM-Newton 
observations of XTE J1810-197 listed in Tabled 

2. XMM-Newton OBSERVATIONS AND RESULTS 

The first XMM-Newton Target of Opportunity (TOO) 
observation of XTE J1810-197 was obtained on 2003 
September 8 and described in Paper I. A second XMM- 
Newton TOO observation of XTE J1810-197 was per- 
formed on 2004 March 11. We use the data obtained 
with the Europea n Photon Imaging Camera (EPIC, 
iTurner et all 12003). EPIC consists of three CCD cam- 
eras, the EPIC pn and the two EPIC MOS imagers, 
which are sensitive to X-rays in 0.1 — 10 keV nominal 
range. The observing modes were different for the two 
observations. In 2003 September, one of the EPIC MOS 
cameras was operated "small window" mode, for which 
the field of view (FOV) of the central CCD was reduced 
to 1.'8 x 1'.8 and read out in a time of 0.3 s. The second 
EPIC MOS camera was operated in "full frame" mode, 
which is read out in 2.7 s integrations over the 30' FOV. 
The EPIC pn was operated in "small window" mode, 
which provides 6 ms time resolution over a 4.'3 x 4f3 FOV. 
In 2004 March, both of the EPIC MOS cameras were op- 
erated in "small window" mode, while the EPIC pn was 
operated in "large window" mode, which provides 48 ms 
time resolution over a 13.'5 x 26' FOV. Thus, all data sets 
except the "full frame" MOS in 2003 September were of 
high enough time resolution to study the pulse profiles 
of XTE J1810-197. In addition, XTE J1810-197 was 
detected serendipitously in a short observation, as yet 
unpublish ed, from the X MM-Newton Galactic Plane Sur- 
vey ijHands et alJl2004j) on 2003 October 12. We include 
results from this survey observation for completeness, in- 
cluding timing information from the pn camera in "full 
frame" mode with 73.4 ms time resolution. 

In order to make the cleanest possible comparison 
among the XMM-Newton observations, we extracted all 
calibrated photon event data files using the latest ver- 
sion of the XMM-Newton Science Analysis System [re- 



lease version SAS 6.0.0 (20040309.1146)], by applying 
the processing chains to each set of Observation Data 
Files (ODF). Screened photon event lists were generated 
by applying the standard XMM-Newton good-time filter 
criteria and selecting CCD PATTERN < 4 for our EPIC 
pn spectral analysis. A total of 11.5 ks, 6.9 ks, and 17.0 ks 
of good EPIC pn exposure time were acquired during the 
three observing epochs, which translates to an effective 
live time of 8.1 ks, 6.2 ks, and 15.8 ks. Resulting back- 
ground subtracted count rates were 10.6 s _1 , 4.8 s _1 , and 
5.8 s _1 , where the second one is affected by vignetting at 
the off-axis angle of 8.'8. The observations were mostly 
uncontaminated by flare events. Photon arrival times 
were converted to the solar system barycenter using the 
Chandra derived coordinates of the source given in Pa- 
per I, (J2000) 18 h 09 m 51. s 08, -19°43'51."7. We manually 
corrected a jump of +1 s in the EPIC pn photon arrival 
times that affected part of the 2004 March observation. 
Timing anomalies are frequent in the pn data; evidently, 
they have not been completely eliminated in SAS 6.0.0. 

2.1. Spectral Modeling 

In Table [21 we summarize the results of spectral 
analysis of previous and new observations. Although 
of relatively poor quality, the quiescent ROSAT spec- 
trum of XTE J1810-197 can be reasonably well fitted 
with a blackbody of kT = 0.18 ± 0.02 keV covering 
1.2 x 10 13 cm 2 , or 2/3 the area of a 12 km radius NS 
at d = 5 kpc. In contrast, the X-ray spectra in outburst, 
obtained from the XMM-Newton TOO observations, re- 
quired two components, which is typical of AXPs (e.g., 
iMarsden fc White! l200l|) . We fitted XMM-Newton data 
from the EPIC pn CCD only. The fast read out of this 
instrument ensures that its spectrum is not affected by 
photon pile-up. Source spectra were accumulated in a 
45" radius aperture which encloses > 95% of the encir- 
cled energy. Background was taken from a circle of the 
same size displaced 2 '3 along the readout direction. The 
spectra were grouped into bins containing a minimum 
of 400 counts (including background) and fitted using 
the XSPEC package. The "wabs" photoelectric absorp- 
tion cross sections of iMorrison &: McCammonl |l983), 
not including Thomson scattering, were used with the 
lAnders fc Ebiharal l)1982f ) abundances. For each model, 
the column density was treated as a single parameter in 
a linked fit to all three epochs, after verifying that inde- 
pendent fits yielded consistent values of TVjj. The other 
fitted parameters are also consistent with their values 
obtained in fits with independent Ah. 

The XMM-Newton spectra are equally well fitted by 
a power-law plus blackbody model, as commonly quoted 
for AXPs, or by a two-temperature blackbody model that 
we will argue is more physically motivated. Sample fits 
are shown in Figure |2 But the ambiguity about which 
X-ray spectral model is appropriate affects our estimate 
of the distance to the source by contributing systematic 
uncertainty to the fitted column density. The X-ray mea- 
sured Nn = (1.02 ± 0.05) x 10 22 cm~ 2 in the case of the 
power-law plus blackbody model is significantly differ- 
ent from the value of (0.65 ± 0.04) x 10 22 cm~ 2 in the 
double blackbody model (see Table|2Jl. We argued in Pa- 
per I that, in cither case, XTE J1810— 197 is likely to be 
closer than the neighboring supernova remnant Gil. 2- 
0.3, which has an H I absorption kinematic distance of 
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TABLE 1 

XMM-Newton Observations of XTE J 1810- 197 





EPIC pn 


EPIC pn 


EPIC pn 


Parameter 


2003 Sep 8 


2003 Oct 12 


2004 Mar 11 


Off-axis angle (arcmin) 


0.0 


8.8 


0.0 


Exposure time (ks) 


11.5 


6.9 


17.0 


Live Time (ks) 


8.1 


6.2 


15.8 


Count Rate (s _1 ) a 


10.6 


4.8 


5.8 


Epoch (MJD/TDB) b 


52890.5642044 


52924.0000320 


53075.4999960 


Period (s) c 


5.539344(19) c 


5.53943(8) 


5.539425(16) c 



a Background subtracted EPIC pn count rate corrected for detector dead- 
time. 

b Epoch of phase zero in Figure l4l 

c Includes EPIC MOS data. 95% confidence uncertainty in parentheses. 



TABLE 2 
Spectral Fits and Fluxes 



Model 


ROSAT /PSPC 


XMM-Newton/EPIC pn 


Parameter 


1992 Mar 7 2003 Sep 8 


2003 Oct 12 2004 Mar 11 


Blackbody 


Nu (10^ cm" 2 ) 


0.63 (fixed) . . . 




kT BB (keV) 


0.18 ±0.02 




BB Area (cm 2 ) 


1.2 x 10 13 . . . 




Flux (ergs cm -2 s _1 ) a 


6.9 x 10" 13 . . . 




L BB (bol) (ergs s" 1 ) 


1.3 x 10 34 . . . 




Xl(dof) 


1.1(13) 




Power-law + Blackbody 



N H (10 22 cm~ 2 ) b 

r 

kT BB (keV) 
BB Area (cm 2 ) 
PL Flux (ergs cm -2 
Lph (ergs s" 1 ) 
BB Flux (ergs cm -2 
L BB (bol) (ergs s" 1 ) 
Total Flux (ergs cm" 
xg(dof) b 



-i )a 



1.02 ±0.05 
3.87 ±0.25 
0.67 ±0.01 

6.2 x 10 11 
1.17 x 10" 11 

2.7 x 10 35 
2.80 x 10" 11 

1.3 x 10 35 
3.97 x 10- 11 

1.077(187) 



1.02 ±0.05 
3.51 ±0.15 
0.68 ±0.01 
5.2 x 10 11 

1.42 x 10" 11 
3.6 x 10 35 

2.47 x lO" 11 
1.1 x 10 35 

3.89 x 10- 11 
1.077(84) 



1.02 ±0.05 
3.75 ±0.20 
0.68 ±0.01 

2.7 x 10 11 
8.47 x 10~ 12 

1.8 x 10 35 
1.31 x 10" 11 

6.0 x 10 34 
2.16 x lO" 11 
1.077(194) 



Double Blackbody 
N n (10 22 cm~ 2 ) b 
fcTi (keV) 
kT 2 (keV) 
BB1 Area (cm 2 ) 
BB2 Area (cm 2 ) 
BB1 Flux (ergs cm" 2 s 
BB2 Flux (ergs cm" 2 s 
iBBi(bol) (ergs s" 1 ) 
L B B2(bol) (ergs s _1 ) 
Total Flux (ergs cm -2 
X^(dof) b 



-1)3 
-1)3 



-1)3 



0.65 ±0.04 
0.26 ± 0.02 
0.68 ±0.02 

1.1 x 10 13 
6.4 x 10 11 

4.24 x 10" 12 
3.51 x 10" 11 

5.2 x 10 34 
1.4 x 10 35 

3.93 x 10" 11 
1.076(187) 



0.65 ±0.04 
0.29 ± 0.04 
0.71 ±0.03 
6.6 x 10 12 
5.1 x 10 11 
5.37 x 10~ 12 
3.03 x 10" 11 
5.1 x 10 34 
1.3 x 10 35 
3.84 x 10~ n 
1.076(84) 



0.65 ±0.04 
0.27 ± 0.02 
0.70 ±0.01 

6.8 x 10 12 

2.9 x 10 11 
3.52 x 10~ 12 
1.78 x 10" 11 

3.9 x 10 34 
7.2 x 10 34 
2.13 x 10" 11 
1.076(194) 



Note. — Uncertainties are 90% confidence for two interesting parameters. 
a Absorbed flux in the 0.5-10 keV band. 
b Parameter derived from a linked fit to all epochs. 
c Unabsorbed luminosity in the 0.5—10 keV band assuming d = 5 kpc. 



5 kpc ijGreen et al.1 119881) and an X-ray measured co l- 
umn density of « 1.4 X 10 22 cm -2 l|Vasisht et al.lll996l) . 
With the distance somewhat uncertain, we parameterize 
derived quantities in terms of c/5, an upper limit on the 
distance, in units of 5 kpc. 

Considering first the 2003 September observation, the 
best fit to the power-law plus blackbody model has pho- 
ton index T = 3.87 ± 0.25 and kT BB = 0.67 ± 0.01 keV, 
with a fit statistic of xl — 1-077 for 187 degrees of free- 



dom. The flux for each of these components is given in 
Table |21 The blackbody flux is equivalent to a spherical 
area of 6.2 x lO 11 ^ 2 cm 2 , or w 3% of the surface area 
of a NS. It is noteworthy that the 2004 March spectrum 
had nearly identical parameters, with T = 3.75±0.20 and 
kTsB — 0.68±0.01 keV, while the fluxes of the power law 
and blackbody declined by 28% and 53%, respectively. 

A two-tcmperature blackbody model yielded a fit 
of essentially the same quality, with a fit statistic of 
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Fig. 2. — XMM-Newton EPIC pn spectra of XTE J1810— 197 from the earliest and latest epochs, fitted with a power-law plus blackbody 
model or a double blackbody model as described in the text and in Table |21 Also shown are the residuals from the best-fit models, (a) 
Power-law plus blackbody fit in 2003 September. (6) Double blackbody fit in 2003 September, (c) Power-law plus blackbody fit in 2004 
March, (d) Double blackbody fit in 2004 March. 



\l = 1.076 for 187 degrees of freedom. The 2003 
September spectrum has kT% — 0.26 ± 0.02 keV and 
kT2 = 0.68±0.02 keV associated with blackbody areas of 
1.1 x 10 13 d\ cm 2 and 6.4 x 10 11 d 2 cm 2 for the cooler and 
hotter temperatures, respectively, which are 60% and 
« 4% of the NS surface. Thus, the power-law component 
of the previous model is replaced with a cooler blackbody, 
while the hotter thermal component retains the same pa- 
rameters in both models. In the double blackbody model, 
the cooler component covers a large fraction of the NS 
surface, but at a temperature significantly higher than 
the quiescent ROSAT measured value. In this parame- 
terization, the cooler and hotter blackbody luminosities 
declined by 25% and 50%, respectively, between 2003 
September and 2004 March. Since their temperatures 
remained the same to within 3% (feTi = 0.27 ± 0.02 keV 
and kT 2 = 0.70 ± 0.01 keV in 2004 March), their decline 
in luminosity is attributable to a decrease in area. 

In order to extract further information about the dis- 
tribution of emission over the stellar surface, we analyzed 
phase-resolved spectra by fitting the two-temperature 
blackbody model to the 2003 September data in 10 phase 
bins. In particular, we kept the temperatures and column 
density fixed at the values listed in Table [5] while fitting 
for the normalization (intensity) in each phase bin. This 
is equivalent to determining the relative projected area 



of emission as a function of rotation phase. This proce- 
dure yields acceptable fits in all phase bins, with reduced 
X 2 values ranging from 0.73 to 1.21. In nine out of 10 
bins, the reduced \ 2 is less than 1.08. Since these fits are 
statistically acceptable as a set, no significant test can be 
made for variations in additional parameters such as the 
temperatures or column density. 

The results of the phase-resolved spectral fitting are 
shown in Figure [3] As can be expected, the intensity of 
the hot component varies by a large factor, « 3.5, while 
the larger cool component varies by only a factor of 1.7. 
Both components peak at the same phase. This is con- 
sistent with the energy-dependent light curves (Fig. QJ, 
in which the pulsed fraction increases with increasing 
energy. The relative amplitudes and phases of the two 
components are therefore consistent with the picture of 
a small, hot region surrounded by a cooler, concentric 
annulus that occupies of order half the surface area of 
the star. Neither component goes to zero at any rotation 
phase. In particular, if the hot component were com- 
pletely eclipsed, the spectral decomposition in Figure [5] 
indicates that the light curves at E > 3 keV should go 
to zero, which clearly they do not. 

While the overall fits of the two models are of equal 
quality, the double blackbody model does leave signif- 
icant residuals above 6 keV, which may be indicative 
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Fig. 3. — Comparison of the two-component blackbody intensi- 
ties as a function of rotation phase in the 2003 September XMM- 
Newton EPIC pn spectrum. The temperatures kTi and KT2 are 
held fixed at the values listed in Table |2] while the normalization 
constants are fitted. Inset: The ratio of blackbody normalization 
constants as a function of rotation phase. 



of an additional harder component, especially in the 
2004 March spectrum. We do not have enough high- 
energy information to model this excess quantitatively, 
but it may represent nonthermal magnetospheric emis- 
sion similar to the hard X-ray co mponent recently re- 
ported from the AXPs IE 1841-045 llBassani et aU2004 
iKuiper. Hermsen. fc Mended 120041: iMolkov et all 12004') , 
1RXS J17084 9.0-400910 lIRevnjvstey et al .1 120041) "^^ 
4U 0142+614 l|den Hartog et alJl2004|) in RXTE and IN- 
TEGRAL data. 

2.2. Pulsed Light Curves and Period Derivative 

The quiescent ROSAT source was not modulated, with 
a pulsed fraction upper limit of 24% (Paper I). Since 
the first detection of pulsed emission in early 2003, the 
pulsed fraction remained fairly constant at w 54% as 
seen in the two Chandra observations on 2003 August 27 
and November 1 (Paper I). The folded light curves from 
the 2003 September 8 XMM-Newton observation (Fig.EJ) 
show that the pulse peak is somewhat narrower than a si- 
nusoid, and the pulsed fractions increase smoothly with 
energy from 38% at less than 1 keV to » 55% above 
5 keV. As discussed in §2.1, the fact that the high-energy 
light curves do not dip to zero implies that the hot black- 
body component is never completely out of view. In 2003 
October and 2004 March, the XMM-Newton light curves 
have not changed much. The pulse shape and pulsed 
fractions are largely unchanged. The marginal appear- 
ance of narrow components at the peak of the pulse pro- 
file in the highest energy bin may be additional evidence 
of a hard nonthermal contribution. 

The measured barycentric pulse period was 5. 539344 ± 
0.000019 s on 2003 September 8, and 5.539425 ± 
0.000016 s on 2004 March 11 (Table [TJ, from the com- 
bined EPIC pn and MOS data. The measurement from 
the shorter pn observation of 2003 October 12 has a much 
larger error bar, but is consistent with the other values. 
The errors are the 95% confidence level determined from 
the Z\ test. The two precise period measurements can 
be compared to yield P = (5.1 ± 1.6) x 10~ 12 s s _1 over 
the 6 month interval, which implies a characteristic age 
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Fig. 4.— Energy-dependent pulse profiles of XTE J1810-197. 
(a) EPIC pn on 2003 September 8. (ft) EPIC pn on 2004 March 11. 
The epochs of phase zero are given in Table HI Pulsed fractions, 
defined as the fraction of the counts above the minimum in the 
light curve, are indicated in each panel. Background has been 
subtracted. Pulse profiles from 2003 October 12 (not shown) are 
similar to these. 



t c « 17, 000 yr, surface magnetic field B s » 1.7 X 10 14 G, 
and spin-down power E w 1.2 x 10 33 er g s~\ compa- 
rable to the earlier values measured by llbrahim et alJ 
However, it is evident that P is not stable, 



11200. 

having varied by at least a factor of 2 on timescales of 
months, llbrahim et alJ l|2004[) fitted values of P in the 
range (1.1 — 2.2) x 10 -11 s s -1 in the first 9 months of the 
outburst. Increases in P o f order unity have been seen in 
the AXPs IE 2259+586 l|Iwasawa. Kovama. fc HalperrJ 
119921 iKaspi et alJ l200l and 1RXS J170849.0-400910 
Ipall'Osso et all l2003V following their glitches. Even 
larger changes in P ar e seen in the SGRs 190 0+14 
lIKouveliotou et alJll99& IWoods et alJlilM 12001 and 
1806-20 ijWoods et al.ll2000Ll2002|l . although these were 
apparently not associated with episodes of bursting. It is 
possible that XTE J1810— 197 experienced a glitch and 
an increase in P at the time of its outburst, and that P 
is now relaxing to its long-term value. Alternatively, the 
spin-down could have been enhanced in the early stages 
of the outburst by a particle wind that is now declining. 
Sufficiently energetic Alfven waves drive winds that can 
increase the magnetic field strength at the light cylinder 
over the dipole value, thereby inc reasing the spin-down 
torque. l|Thompson fc Blaeslll99S|) . 

We note here that two discrepant period measure- 
ments have been reported from independe n t analy ses of 
the 2004 March observation. iRea et all (l2004af) gave 
P = 5.53974 + 0.00005 s, and later IRea et all <|2004bh 
gave P = 5.539917 ± 0.000005 s. Both of these are 
evidently affected by the error in the EPIC pn tim- 
ing that we discovered in the data (§2). Consequently, 
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the much larger P reported by iRea et all (2004b) is 
not valid. Our correct period measurement in Table ^ 
P = 5.539425 ± 0.000016 s, is confirmed by establish- 
ing consistency, both in period and in phase, among the 
three EPIC instruments' data analyzed separately. 

2.3. Long-term Flux Decay 

llbrahim et alJ l)2004j) fitted the decline of the flux from 
XTE J1810-197 as observed by RXTE with either an 
exponential of time constant 269 ± 25 days, or a power 
law of r 13 where 0.45 < (3 < 0.73. The allowed range of 
[3 comes from the uncertainty in the initial epoch of the 
outburst, 2003 November 17 < t < 2004 January 23. 
The new XMM-Newton observations extend the baseline 
by an additional 5.5 months, although it is difficult to 
make a precise comparison between RXTE and XMM- 
Newton fluxes because of the very soft spectrum of the 
source, and the considerable background and source con- 
fusion that plagues RXTE in this region of the Galac- 
tic plane. Nevertheless, the first XMM-Newton obser- 
vation is consistent in flux w ith the contempora neous 
RXTE monitoring (Fig. HI and llbrahim et aLll2004^ . and 
the September-March XMM-Newton points alone fit an 
exponential time constant of 300 days, marginally consis- 
tent with the RXTE determined decay. However, if we 
fit the XMM-Newton points with a power-law decay, then 
the allowed exponent range is 1.03 < (3 < 1.26, quite a 
bit steeper than the RXTE fitted rate. Continued mon- 
itoring is thus able to discriminate among exponential, 
power law, or more complicated decays. It already seems 
that an exponential is a better fit than a power law. 

The total energy emitted in the outburst can be esti- 
mated by integrating the fitted exponential decay, i.e., 
multiplying the peak luminosity by the decay constant. 
Assuming a peak bolometric luminosity of 6 x 10 35 <i§ 
ergs s" 1 at the onset of the RXTE observered outburst 
ijlbrahim et alJ 120041) and a time constant of 300 days 
yields a fluence of 1.5 x 10 43 d\ ergs. 

3. DISCUSSION 

3.1. Problems with a Power-Law X-ray Spectrum 

A two-component X-ray spectral fit is commonly re- 
quired for AXPs, but its physical significance is poorly 
understood. While it is conventional to fit a power- 
law plus blackbody model, the interpretation of the re- 
sulting steep power law, e.g., T w 3.8 in the case of 
XTE J1810— 197, is not clear. When employing power 
laws as a means of fitting excess high-energy emission, 
problems arising from the dominance of such power laws 
at low energy are usually ignored. Here, we point out 
that such a steep power law cannot be c o nnect ed to the 
faint IR fluxes measured bv llsrael et alJ l)2004j) without 
assuming severe and perhaps unphysical cut-off mecha- 
nisms. 

Accor ding to Figu r e EK, which is equivalent to Fig- 
ure 2 of llsrael et alJ ((2004) , the if-band measurement 
and 0.5 keV flux density of XTE J1810-197 as predicted 
from the power-law fitted component can be bridged by 
a power law of F v oc v a , where a = +0.6. Notably, 
this value of a is in excess of the low-frequency limit 
of an optically thin synchrotron spectrum, which has 
F v oc 2/+ 1 ' 3 . Alternatively, the standard synchrotron 
self-absorption mechanism, which allows a low-frequency 



spectrum F v oc j/ +5 / 2 , can be invoked to avoid exceeding 
the observed infrared flux. But this would require the 
self-absorption frequency v a to fall in the narrow (and 
unobservable) range 10 16 — 10 17 Hz. It also places a re- 
quirement on the radius r of the synchrotron source, 

r = 6 x 10 14 B 1/4 ^ 5 / 4 Fll 2 d cm, 

where F Va is the flux density at v a , extrapolated from the 
fitted X-ray power law. This relation assumes a uniform 
B field unlike the magnetosphere of a neutron star. But 
if we combine it with an estimate of the local dipole field 
at radius r, 




with B s « 2 x 10 14 G, then the synchrotron source must 

have r < 3 x 10 7 dtj 7 cm. Such small values of r are 
inconsistent with the synchrotron assumption since the 
associated magnetic field B(r) > 10 10 G is large enough 
that the soft X-ray spectrum should instead be a broad 
cyclotron feature. These serious problems with a syn- 
chrotron model cast doubt on the utility of a steep power- 
law fit. 

Alte rnatively, while FrhomDSon . Lvutik ov. fc Ku lkarni 
(2002) propose that multiple resonant Compton scatter- 
ing of thermal photons can generate a high-energy tail, 
this mechanism does not address the origin of the re- 
quired low-energy excess, which is just as much a part of 
the steep power-law component as the high-energy tail. 
Therefore, we conclude that a power law is not a correct 
description of the actual spectrum. While it fits ade- 
quately in the limited X-ray band, it has no compelling 
physical explanation and is in conflict with lower-energy 
measurements. Unfortunately, the large column densi- 
ties to most AXPs prevent an accurate characterization 
of their intrinsic soft X-ray spectra because the correc- 
tions for absorption become highly uncertain when they 
are large. 

The blackbody component in the power-law plus black- 
body model is easier to accept since its effective area is 
initially « 6 x 10 11 d\ cm 2 (neglecting unknown geometric 
and beaming factors), decreasing to « 3 x 10 11 d\ cm 2 in 
the third observation. Because these areas are both small 
compared to the surface of the NS, the large pulsed frac- 
tion of this component, which dominates the flux from 
2-5 keV and does not change as the flux declines, can 
be understood as rotational modulation of a surface hot 
spot. 

The pulsed light curves provide an independent diag- 
nostic of emission mechanisms. The strict pulse-phase 
alignment and subtle change of pulse shape as a function 
of energy in Figure^lis unnatural under the hypothesis of 
two entirely different emission mechanisms and locations. 
The power-law component of the spectral fit dominates 
below 1.5 keV while the blackbody component dominates 
between 2 and 6 keV. The observed monotonic increase in 
pulsed fraction as a function of energy is not accounted 
for in such a spectral model. This problem applies to 
AXPs in general, as discussed in detail by Ozel, Psaltis, 
& Kaspi (2001). A hard X-ray tail can arise from radia- 
tive transfer effects in a realistically modeled NS atmo- 
sphere ijOzell 120011 lOzel et al.ll200lt iPerna et alJl200l|) . 
thus reducing the need for a power-law component. This 
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Fig. 5. — Broad-band spectrum of XTE J1810— 1 97 as derived from fits to the XMM-Newton EPIC pn data of 2003 September, and 
contemporaneous infrared and optical data points of llsrael et al.l 120041) from 2003 October. Filled circles are observed points, and open 
circles are corrected for extinction using the X-ray fitted Nh from each model and the Pre dehl fc Sch mitt . (.1995) relation Na/Ay = 1.8x 10 21 
cm -2 mag -1 . Fitting parameters of the two-component X-ray spectra are taken from Tabled] The unfolded, unabsorbed spectrum is 
represented here, (a) Power-law plus blackbody fit to the X-ray spectrum, with Nn = 1.02 X 10 cm -2 and Ay = 5.7 mag. This panel is 
equivalent to Fig ure 2 of llsrael etUl 120041) . (ft) Double blackbody fit to the X -ray spectrum, with TVjj = 0.65 X 10 cm and Ay = 3.6 
mag. 



raises the possibility that a purely thermal model may 
account for the spectra and pulse profiles. 

3.2. Advantages of a Purely Thermal X-ray Spectrum 

In the alternative double blackbody spectral model, a 
large area of ~ 60% the surface of the NS is implied for 
the cooler blackbody component, while the hotter com- 
ponent has an initial area of w 4% of the NS surface. 
In its favor, a purely thermal model can explain why 
the pulsed fraction increases with energy, and why the 
pulses are in phase at all energies, assuming the geome- 
try to be that of a small hot spot surrounded by a cooler 
annulus. The phase-resolved spectroscopy discussed in 
§2.1 supports this picture. Also, the extrapolated X-ray 
spectrum in this model does not exceed the observed IR 
fluxes (Fig. Eb). IR and optical emission from AXPs, 
including XTE J1810— 197, are too bright to be surface 
thermal emission, but are likely to be magnetospheric 
synchrotron radiation. Unlike the X-rays, the optical/IR 
emission does not exceed the spin-down luminosity, and 
could be rotation-powered in quiescence (jOze]H2D0^l . al- 
though additional energy deposited by magnetic field de- 
cay could increase its brightness during outbursts. 

Parenthetically, we mention the case of the AXP 1RXS 
J170849.0-400910. A double blackbody was also one 
of the models fitted to the X-ra y spe ctrum of 1RXS 
J170849.0-400910 bv llsrael et all l|2001ft . with an inter- 
pretation similar to ours. On the other hand, 1RXS 
J170849. 0-400910 is unique among AXPs for having a 
strong energy dependence of its pulse phase, shif t ing by 
w 0.3 cycles between 2 and 10 keV. iRea et alJ (2003) 
interpret this as most likely due to a phase dependence 
of a power-law spectral index. Interestingly, this is the 



only AXP for which good evidence of a cyclotron ab sorp- 
tion feature is found, at 8.1 keV l)Rea et al J 12003ft . and 
therefore a true deviation from blackbody emission. 

While a purely thermal model may be favored by the 
simple arguments presented above, detailed modeling of 
the light curves will be needed to verify that the soft- 
est X-rays from XTE J1810— 197, which come from the 
cooler, larger area blackbody, can have a pulsed fraction 
as high as the observed 38%. The emitting area might be 
even larger if realistic atmosphere models are fitted in- 
stead of pure blackbodies ijPerna et~alll200l . One pos- 
sibility is that d is considerably less than 5 kpc, thus 
reducing the emitting area. Another important consid- 
eration is the effect of realistic radiative transfer in highly 
magnetized atmospheres. IQzel et ahl l|2001ft showed that 
the sorts of pulsed light curves that are observed from 
XTE J1810-197 and other AXPs are inconsistent with 
two antipodal hot spots on the surface. However, they 
showed that a single hot spot is capable of producing 
large pulsed fractions as well as the observed energy de- 
pendence from a nearly orthogonal rotator if the mag- 
netic field is assumed to be normal to the surface. [Note 
that lOzel et"aTI pOOl) use a different definition of pulsed 
fraction from the one that we use in Figure^]. The large 
predicted pulsed fractions rely on the "beaming" effect 
of anisotropic radiative transfer in a large magnetic field, 
in which the electron scattering cross sections of both 
the ordinary and extraordinary modes of polarization are 
minimized where the angle between the wave vector k 
and the magnetic field B is zero. 

The quiescent ROSAT spectrum of XTE J1810-197 
is fitted by a simple blackbody comparable to the full 
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surface area of a NS. Note that the temperature mea- 
sured by ROSAT using a single blackbody model is even 
lower that found for the cooler component of the two- 
temperature XMM-Newton fit. This perhaps accounts 
for the failure to detect rotational modulation in the qui- 
escent state. Thus, it is possible that the X-ray mani- 
festation of the outburst of a TAXP is largely due to 
a sudden thermal heating event in the crust, while the 
long-term quiescent flux is powered by slow dissipatio n of 
the core magnetic field ijThompson fe Duncanl [T996I) or 
simply by residual cooling in the pres ence of an insulat- 
ing at mosphere of hydrogen or helium <)Hevl fc HernauisH 

E223>. 

3.3. Afterglow Following Deep Crustal Heating 

Similar to XTE J1810-197, a decay of the thermal 
component in the AXP IE 225 9+586 was seen dur ing the 
bursting episode in 2002 June ijWoods et alJl2004^ . when 
a small, hot region of kT-QB ~ 1.7 keV and R w 1 km 
cooled within one day, leaving a quiescent thermal com- 
ponent of fcTBB ~ 0.5 keV and R rj 6 km. After 1 
day, the 2-10 keV X-ray flux decayed further as a slow 
£-0.22 p 0wer i aw f or a t least 200 days. Hot spots due to 
deep crustal heating by localized magnetic field decay, or 
surface bombardment by accelerated particles, can plau- 
sibly produce the observed spectra and pulse profiles. 
However, heating by magnetic field decay is preferred 
over particles accelerated only on open field lines in the 
dipole geometry, since the canonical dipole polar cap area 
is only 2 x 10 s (P/5.5 s) _1 cm 2 for such slow rotators as 
t hese AXPs, much sm aller than the detected areas. 

Eich ler et al.l l)2003|) reviewed mechanisms for "after- 
glows" from cooling of the crust following SGR like out- 
bursts in magnetars. Deep crustal heating that is sub- 
ject to a long cond uction time to th e surfa ce can plausi- 
bly be the cause. iThomoson et al.l l|2002j) showed that 
the resulting cooling would follow a t~ - 7 power law. 
The initial heating could take place within 10 4 s due 
to large-scale shearing of the crus t by m agnetic stresses. 
iLvubarskv. Eichler. fe Thompson! l)2002|) assumed a de- 
position of ~ 10 25 ergs cm -3 to a depth of 500 m, which 
is < 1% of the magnetic energy density at the surface. If 
occurring over the entire surface, this is ~ 1 x 10 43 ergs, 
comparable to the total X-ray energy emitted during the 
decay of XTE J1810— 197. If such releases occur with re- 
currence time > 20 yr, AXPs could power transient out- 
bursts for > 5000 yr. However, this simple rate estimate 
does not t ake into account heat co nduction into the core 
of the NS. ILvubarskv eTaD IpJOlj) find that only 20% of 
the deposited heat is radiated during the afterglow; the 
rest is conducted into the interior and contributes to the 
average quiescent X-ray or neutrino emission. The ob- 
served quiescent luminosity of 1.3 x 10 34 ergs s _1 could 
be powered by ~ 7 x 10 45 ergs of magnetic energy dissi- 
pated in the core over the presumed ~ 17, 000 yr age of 
the NS. 

This theory has been applied to the slow cool- 
ing observed from SGR 1900+14 lllbrahim et alJ 120011: 
IWoodset aimOOlULenters et al.ll2003l) and SGR 1 627- 
41 (iKouveliotou et all 120031 ILvubarskv et al.ll2002j) fol- 
lowing their outbursts. In the spectra of SGR 1900+14, 
a blackbody component was identified that has a high 
initial temperature, w 2 keV, and radius of w 1.7 km, 
much smaller than a neutron star, showing that heating 



(or heat conduction or radiation) can indeed be local- 
ized. Therefore, it seems that the long timescale decay 
of the flux from XTE J1810-197 could be due to the 
same process that occurs in the outbursts of SGRs. 

While there are still too few measurements to char- 
acterize reliably the long-term decay of the outburst of 
XTE J1810— 197, it does seem already to be declining 
faster than the t~ - 7 power law predicted in the crustal 
heating model. We find a decay index of 1.03-1.26 in the 
interval 9-15 months after the presumed outburst, but 
an exponential of time constant 300 days is a better de- 
scription of the entire light curve since outburst. Similar 
deviations of the pos t-outburst decay curve of SGR 1627- 
41 from a power law ijKouveliotou et all200 3') complicate 
the crustal heating interpretation. Of course, the existing 
cooling models are based on solving a onc-dimcnsional 
heat transfer equation, and do not involve localized heat- 
ing and possible changes in emitting area such as are 
suggested by the observations of XTE J1810-197. 

Another variable AXP for which detailed spectral anal- 
ysis at different intensity levels has bee n performed is 
IE 1048.1-5937. iMereghetti eTaD l|2004|) found that its 
spectral shape remained invariant between states of lu- 
minosity that differed by a factor of 5. In a fitted two- 
component model, the power-law index as well as the 
blackbody temperature remained the same, which im- 
plies that it is not realistic to assume two physically 
distinct emission mechanisms. Instead, the variability 
in flux could be explained in a purely thermal model 
by a change in emitting area. Such an interpretation 
would be consistent with the pulsed fraction of the light 
curve of IE 1048.1-5937, which decreased as the flux in- 
creased ijMereghetti et al.ll2004ft . We interpret the decay 
of XTE J1810-197 along the same lines. 

4. CONCLUSIONS AND FUTURE WORK 

We used three XMM-Newton observations of 
XTE J1810-197 obtained during its decline from 
outburst to make a number of quantitative measure- 
ments, although their robustness needs to be tested 
with additional follow-up observations. The decay is 
consistent with an exponential of time constant « 300 
days. The pulsar spin-down rate is quite variable, 
having decreased to (5.1 ± 1.6) x 10~ 12 s s _1 after 
initial measurements by RXTE found it to vary in 
the range (1 — 2) x 10~ n s s" 1 . The inferred char- 
acteristic age r c w 17, 000 yr, magnetic field strength 
B s « 1.7 x 10 14 G, X-ray spectra, and pulsed light 
curves, are all consistent with the behavior of an 
anomalous X-ray pulsar, albeit one that has varied in 
flux by two orders of magnitude. Assuming that the 
decay continues at the observed exponential rate, the 
total energy radiated in the present outburst will be 
w 1.5 X 10 43 d\ ergs. This is comparable to the heat that 
can be deposited in the crust by conversion of < 1% 
of the magnetic energy density at the surface. Such a 
magnetar mechanism could power outbursts similar to 
the current one at > 20 yr intervals for a large fraction of 
the characteristic age of the pulsar, while the quiescent 
luminosity observed in a number of archival X-ray 
observations could come from the long-term average 
magnetic decay heating of the NS core. Models of deep 
crustal heating typically predict t~ ' 7 power-law decays, 
which are not matched by the sparse measurements 
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so far obtained from XTE J1810— 197. However, more 
complex decay behavior has been seen in several other 
magnetars, which may require refinements to the initial 
simple models, in particular, generalization to two 
dimensions. 

The X-ray spectrum and pulse profiles of 
XTE J1810-197 are typical of AXPs. We discussed the 
problems with allowing a power-law plus blackbody fit, 
which results in a steep power law that is difficult to 
reconcile with the faintness of the optical/IR counter- 
part. The lack of variation in pulse profile as a function 
of energy also argues against the action of two differ- 
ent emission mechanisms. Instead, a purely thermal 
model for the spectrum, here adequately modeled as 
two temperatures, is more physically motivated. The 
required blackbody areas are less than or equal to the 
surface area of a neutron star, and the pulse modulation 
can probably be understood in terms of anisotropic 
heat conduction and radiative transfer in the strong 
magnetic field. Two dimensional cooling models will be 
required to explain why the decay is characterized as a 
decline in emitting area rather than temperature, the 
latter remaining constant, at least initially. The spectral 
data obtained so far are limited in that the existing 
XMM-Newton observations span slightly less than a 
factor of 2 in luminosity. The most stringent constraints 
on models may yet come from X-ray spectra and pulse 
profiles to be obtained over the coming years, as the 
luminosity of XTE J1810-197 is still a factor of 20 - 30 
above its historic quiescent level. 

The existence of serendipitous archival data allowed 
us to determine that XTE J1810— 197 lacked pulsations 
in the quiescent state preceding the current outburst. 
This suggests that CCOs may be interpreted as qui- 
escent AXPs. No periodic signals have been detected 
so far from the CCOs in the Cas A, Puppis A, and 
Kes 79 SNRs, for example. Cas A is known to be a 
very young, w 300 yr object; its c entral source has a 
spectrum consistent with an AXP l|Chakrabartv et al.l 
l200lUMereghetti. Tiengo. fc Israe]ll2002fl . A similar ob- 
ject may be the po int source at the center of Kes 79 
ijSeward et alJ l2003|) . The radio quiet neutron star 
RX J0822-4300 in Puppis A also has a spectrum sim- 



ilarto that of XTE J1810-197 (iBecker fc Aschenbachl 
12002ft . iCaensler. Bock, fc Stapnersl l|2000fl argue that 
RX J0822-4300 may have AXP-like spin parame- 
ters. The CCO in SNR RCW 103, a non-pulsating 
AXP-like object, has been monitored for over a 
decade and it is found to display marked variabil- 
ity on months to hours time scales. [While a pos- 
sible 6 hr periodic variat ion suggests that this ob- 
ject is a low-mass binary l)Becker fc Aschenbachl 120021: 
iSanwal et al.l l2002j) . the evidence is not conclusive.] 
SGRs are hypothesized to be older (~ 10 4 yr) man- 
ifestations of the A XPs llCottherf. Vasisht. fc Dotanl 
I1999H ICaensler et all l200lF AXPs and SGRs have 
long been considered related phenomena, reinforced 
by the recent detection of SGR-like bursts from 
the two AXPs IE 1048.1-5937 and IE 225 9+586 
ijGavriil. Kaspi. fc Woodsll^ol iKaspi et al.l l2003). The 
various "classes" of young NSs that differ significantly 
from rotation-powered pulsars are phenomenologically 
related, possibly through an evolutionary progression. 

Although TAXPs are relatively rare, their short active 
duty cycle suggests the existence of a larger population of 
unrecognized young NSs. To this end, XTE J1810-197 
provides a unique window into this population, with prior 
measurements in the quiescent state and detailed obser- 
vations during its active, pulsed phase. It is important 
to monitor closely the fading of XTE J1810— 197 during 
this possibly one-time opportunity to study the processes 
of thermal and non-thermal emission to which magnetars 
convert their energy. Measurement of the quiescent spec- 
trum is especially important to help identify the faint 
population of likely missed NSs. 
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